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Non-monotonous  temperature  dependence  of  the  spectral  maximum  of 
photoluminescence  in  CdS/ZnSe  superlattices 

S.  A.  Tarasenko f , A.  A.  Kiselevf , E.  L.  Ivchenkof , A.  Dinger:;:.  M.  BaldaufA, 

K.  Klingshirn$  and  H.  Kaltt 

t Ioffe  Physico- Technical  Institute,  St  Petersburg,  Russia 

t Institute  for  Applied  Physics,  Karlsruhe  University,  D-76128  Karlsruhe,  Germany 

Abstract.  Macro-  and  micro-photoluminescence  (PL)  spectroscopy  is  applied  to  investigate  exciton 
localization  in  cubic  CdS/ZnSe  type-II  superlattices  (SL)  in  the  temperature  range  between  5 K and 
35  K.  The  non-monotonous  shift  of  the  macro-PL  maximum  with  increasing  temperature  reveals 
the  kinetic  contribution  of  acoustic-phonon-assisted  exciton  hopping  processes.  The  kinetic  theory 
of  localized  excitons  has  been  generalized  from  zero  to  finite  temperatures  allowing  for  exciton 
hopping  from  low-energy  to  high-energy  localized  sites  as  well  as  for  the  upwards-downwards 
correlation.  The  results  of  kinetic  theory  are  compared  with  those  of  computer  simulation  and  the 
experimental  data. 


In  this  work  the  properties  of  localized  excitons  of  several  cubic  CdS/ZnSe  superlattices 
have  been  investigated  using  spatially  integrated  (macro-PL)  and  spatially  resolved  (/i- PL) 
photoluminescence  spectroscopy.  The  band  alignment  of  this  heterostructure  is  of  type-II 
with  the  conduction  and  valence  band  off-sets  of  about  800  meV  and  500  meV,  respectively: 
the  CdS  (ZnSe)  layers  form  the  potential  wells  for  the  electrons  (holes)  [ 1.  We  concentrate 
here  on  the  macro-  and  jt-PL  results  of  a 200  period  (19  A/19  A)  CdS/ZnSe  superlattice 
grown  by  molecular-beam  epitaxy  on  GaAs(001)  | ].  The  PL  spectra  discussed  in  the 
following  were  excited  by  the  35 1 nm  line  of  an  Ar+-ion  laser  and  measured  with  a spatial 
resolution  of  about  1 gm  (/<-PL)  or  about  50  //  m (macro-PL). 

The  macro-PL  spectrum  is  displayed  together  with  the  PLE  spectrum  in  the  inset  in 
Fig.  1.  In  Fig.  2(a)  jt-PL  spectra  are  shown,  which  have  been  recorded  at  two  different 
temperatures.  Besides  a broad  PL  background,  narrow  line  emission  is  observed  which  is 
typical  for  strongly  localized  excitons.  Full  circles  in  Fig.  1 show  the  measured  temperature 
dependence  of  the  relative  shift  for  the  macro-PL  maximum 

AEm=ha>m(T)-hcom(0)-8Eg(T),  (1) 

where  com{T)  is  the  spectral  position  of  the  PL  maximum  at  the  temperature  T,  8Eg{T) 
is  the  temperature  shift  of  individual  narrow  lines  in  the  /z-PL  spectra  which  follows  the 
temperature  dependence  of  the  SL  band  gap.  The  value  of  A£m  obviously  represents  a 
purely  kinetic  contribution  to  the  shift. 

For  numerical  simulation  we  generated  two-dimensional,  random,  uniform  distributions 
of  localization  sites  on  a square  cell  of  sufficiently  large  size  with  the  chosen  energy  density 
of  states,  g{e)go  exp  (— e/eo)-  The  PL  spectrum  was  calculated  as  a sum  of  contributions  of 
individual  localized  sites,  each  one  broadened  by  a Lorentzian.  For  the  simulation  of  macro- 
PL  spectra  we  generated  > 2000  configurations  to  obtain  a sufficiently  smooth  spectrum. 
In  order  to  take  into  account  qualitatively  the  inhomogeneous  broadening,  the  macro-PL 
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Fig.  1.  Temperature  dependence  of  the  shift  of  macro-PL  maximum  defined  by  Eq.  ( I ) and  related  to 
the  low-temperature  PL  Stokes  shift  equal  to  18  meV  foraCdS/ZnSe  (19  A/19  A)  SL.  Full  circles, 
experiment;  solid  line,  computer  simulation.  In  inset:  The  PL  and  PLE  spectra  are  compared 
revealing  a Stokes  shift  of  1 8 meV. 


Fig.  2.  (a)  /x-PL  spectra  (spatial  resolution  1 pm)  measured  from  a CdS/ZnSe  (19  A/19  A)  SL 
at  T — 5 K and  20  K.  (b)  Computer-simulated  p-PL  spectra  for  T = 0 K and  T = 22  K.  The  shift 
SEg(T)  is  ignored. 

spectra,  simulated  for  a fixed  value  of  the  exciton  mobility  edge  Eq,  were  convoluted  with 
a Gaussian  describing  the  distribution  of  Eq.  For  the  simulation  of  /x-PL  spectra  (see 
Fig.  2(b))  we  took  200  subsystems  each  containing  100  localized-exciton  sites  randomly 
distributed  with  equal  probability  within  a square  area  and  with  the  weight  g(s)  in  energy. 
The  exciton  mobility  edge  in  each  subsystem  was  chosen  randomly  in  accordance  with  the 
Gaussian  distribution. 

The  parameters  chosen  to  calculate  the  PL  spectra  are  as  follows:  coqxq  — 3xl04, 
TT(a/2)2go£o  = 0.2  (the  notations  are  the  same  as  in  [ ]).  The  inhomogeneous  broadening 
is  described  by  the  Gaussian  with  a width  3.6eo-  At  zero  temperature  the  PL  maximum 
occurs  at  Eq  — hojm  — 3.76sq.  Since  the  PL  Stokes  shift  equals  18  meV  (see  the  inset  in 
Fig.  1)  we  obtain  for  £o  the  value  « 4.8  meV. 

We  have  generalized  the  kinetic  theory  [ , ] of  localized  excitons  from  zero  to  low  but 
finite  temperatures.  At  T — 0 each  site  is  characterized  by  two  parameters:  the  localization 
energy  s — Eq  — E^  ( E\  is  the  localized-exciton  excitation  energy  E)  and  the  distance 
r to  the  nearest  lower-lying  site.  At  finite  temperature  each  site  is  characterized  by  three 
parameters,  s,  s'  and  r,  where  s'  is  the  localization  energy  of  the  optimum  site  for  the  local 
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hopping  process.  If  the  probability  to  occupy  such  a site  is  denoted  by  f(e,  s',  r),  then  the 
energy  distribution  N(s)  which  determines  the  PL  spectral  intensity  is  given  by 

pOO  pOO 

N(e)  = g(s)  / els’  / drPe(s',  r)f(e , e',  r)  . (2) 

Jo  Jo 

Here  P£(s',  r)  is  the  distribution  of  optimum  neighbors  in  energy  and  space.  We  propose 
a rather  simple  and  effective  form  of  the  approximate  kinetic  equation  for  f(e,  s',  r)  and 
write  it  as  [ ] 

1-  w(s,  s',  r)l  f(s,  s',  r)  — (3) 

.to  J ' 

r°°  r°°  »(£,) 

-/  / ds\dr\1———  P£^{s,r\)w{s\,  s,r\)  f {s\,  s,r\)  = T0  . 

Jo  Jo  g(e) 

The  notations  used  are  as  follows:  to  is  the  recombination  time.  To  is  the  generation 
rate  from  delocalized  states,  w{s,  s',  r)  = w(s,  s',  r)p(s,  s',  r),  w(s,  s',  r)  is  the  exciton 
transfer  rate  for  the  e — ► s'  transition  between  the  sites  O and  O'  separated  by  a distance 
r,  for  e < s',  the  factor  p(s,  s',  r ) is  set  to  1 and,  for  e > s'  it  is  equal  to  the  probability 
that  near  the  sites  O,  O'  there  exists  a localization  site  O"  satisfying  the  following  two 
conditions:  (i)  the  upward  hopping  process  O ->  O'  is  more  probable  than  the  process 
O — ► O"  and  (ii)  the  hopping  O'  — >■  O"  is  more  probable  than  the  downward  process 
O'  — > O.  The  factor  p takes  into  account  that,  at  nonzero  but  low  temperatures,  the  major 
part  of  excitons  which  participate  in  the  upwards  transition  0-^>-0'  return  back  from  the 
site  O'  to  the  site  O and  either  recombine  on  the  site  0 or  relax  from  0 to  other  localization 
sites  different  from  O’ . Note  that  in  [ ] the  difference  of  p from  1 is  ignored. 

In  Fig.  3(a)  we  present  the  relative  shifts  of  the  PL  maximum  as  a function  of  temperature 
calculated  neglecting  inhomogeneous  broadening.  The  dashed  curve  is  calculated  by  using 
numerical  solution  of  Eq.  (3)  for  p = 1.  One  can  observe  a remarkable  discrepancy  of  the 
dotted  curve  from  the  solid  curve  which  is  the  exact  result  of  computer  simulation.  The 
dotted  curve  in  Fig.  3(a)  is  calculated  by  using  the  proposed  approach  to  account  for  the 
upwards-downwards  correlation.  It  shows  a good  agreement  with  the  results  of  computer 
simulation.  The  calculated  non-monotonous  behaviour  of  the  Stokes  shift  with  temperature 
was  reported  by  Zimmermann  et  al.  [ ] and  Baranovskii  et  al.  | ].  Physically,  it  can  be 
explained  taking  into  account  that,  at  T = 0,  the  photoluminescence  is  dominated  by 


Fig.  3.  Relative  shift  of  the  macro-PL  maximum  as  a function  of  temperature  calculated  neglect- 
ing (a)  and  taking  into  account  (b)  inhomogeneous  broadening.  Solid  lines,  computer  simula- 
tion; dotted  lines,  the  proposed  kinetic  theory;  dashed  curve,  the  kinetic  theory  [ ] which  ignores 
upwards-downwards  correlation. 
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excitons  finding  themselves  on  accidentally  isolated  localization  sites  acting  as  pores.  For 
such  sites  the  lifetime  with  respect  to  hopping  to  the  nearest  lower-energy  neighbor  exceeds 
the  recombination  time,  to.  At  finite  temperatures,  an  exciton  trapped  by  an  effective  pore 
has  a possibility  to  be  shaken  down  as  it  hops  first  to  the  nearest  higher-energy  neighbor  and 
then  to  another  lower-energy  site.  It  is  the  factor  p(e,  s',  r)  that  estimates  the  probability 
for  this  possibility  to  be  realized. 

Figure  3(b)  presents  the  temperature-induced  shift  of  the  PL  maximum  after  the  PL 
spectrum  is  convoluted  with  the  Gaussian.  The  deviation  between  the  solid  and  dashed 
curves  can  be  interpreted  in  the  following  way.  If  the  PL  spectrum  calculated  for  the 
homogeneous  mobility  edge  is  symmetrical  with  respect  to  the  maximum  then  the  inhomo- 
geneous broadening  makes  no  effect  on  the  spectral  maximum  position.  However,  if  the 
PL  spectral  maximum  calculated  for  a fixed  £o  is  asymmetrical , the  inhomogeneous  broad- 
ening results  in  a shift  of  the  maximum  position.  At  zero  temperature  the  both  calculated 
PL  spectra  almost  coincide  in  shape.  Hence  the  difference  between  the  solid  and  dashed 
curves  in  Fig.  3(b)  appears  because  the  homogeneous  macro-PL  spectral  shape  is  temper- 
ature dependent  and,  at  final  temperatures,  asymmetries  of  the  simulated  PL  spectrum  and 
of  that  calculated  by  using  the  kinetic  theory  are  different. 
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